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ABSTRACT 

OO ■ 

Aims. We report upper limits to the very high energy flux {E > 100 GeV) of the flat spectrum radio quasar 3C454.3 (z = 0.859) derived by 
IL" the Cherenkov telescope MAGIC during the high states of July/ August and November/December 2007. We compare the upper limits derived in 
, both time slots with the available quasi-simultaneous MeV-GeV data from the AGILE y-ray satellite and interpret the observational results in the 
context of leptonic emission models. 

Methods. The source was observed with the MAGIC telescope during the active phases of July-August 2007 and November-December 2007 
and the data were analyzed with the MAGIC standard analysis tools. For the periods around the ends of July and November, characterized by 
■ " " the most complete multifrequency coverage, we constructed the spectral energy distributions using our data together with nearly simultaneous 
multifrequency (optical, UV, X-ray and GeV) data. 

Results. Only upper limits can be derived from the MAGIC data. The upper limits, once corrected for the expected absorption by the extragalactic 
background light, together with nearly simultaneous multifrequency data, allow us to constrain the spectral energy distribution of 3C454.3. The 
data are consistent with the model expectations based on the inverse Compton scattering of the ambient photons from the broad line region by 
relativistic electrons, which robustly predicts a sharp cut-off above 20-30 GeV. 

Key words, quasars: individual (3C454.3) - Gamma rays: observations - Gamma rays: theory 



1. Introduction 

The present list of known extragalactic sources of Very High 
Energy (VHE; defined here as £ > 100 GeV) radiation in- 
cludes 24 source^ (e.g. Aharonian et al. 2008a, De Angelis et al. 
20080 As expected, the majority of these sources (18) belong 

Send offprint requests to: F. Tavecchio, N. Mankuzhiyil and V. Scapin 
' 25, pending the confirmation of the radiogalaxy 3C66B detected by 

MAGIC (Aliu et al. 2008). 
^ see also ht tp : //www . mppmu . mpg . de/~rwagner/sources | 



to the high-peaked BL Lac class. The remaining 6 are four low 
peaked-BL Lac objects (BL Lac: Albert et al. 2007; W Comae: 
Acciari et al. 2008; S5 0716-H71: Teshima et al. 2008; 3C66A: 
Swordy et al. 2008), a radiogalaxy (M87: Aharonian et al. 2003, 
2006a) and a quasar (3C279: Albert et al. 2008a). 

Although the detection of 3C279 indicates that also quasars 
can, to some extent, emit VHE radiation, general theoreti- 
cal arguments support the view that, due to internal absorp- 
tion (e.g. Liu & Bai 2006, Reimer 2007) and/or to the de- 
crease of the cross section for inverse Compton scattering (e.g. 
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Tavecchio & Ghisellini 2008) powerful Flat Spectrum Radio 
Quasars (FSRQs) cannot be prominent VHE emitters. Moreover, 
FSRQs are generally located at relatively high redshifts, imply- 
ing a huge absorption of y-rays by the Extragalactic Background 
Light (EBL). On the other hand, the detection of these sources at 
VHE would be important for our understanding of their structure 
and of acceleration/emission mechanisms and would provide a 
unique opportunity to probe the EBL at relatively high redshifts, 
allowing to study its evolution over cosmic time. 

3C454.3 (z = 0.859) is a well known FSRQ, detected several 
times in the y-ray band by the EGRET telescope onboard CGRO, 
with an average photon index of F = 2.2 (Hartmann et al. 1999). 
In 2005 it underwent a very active phase in optical and X-ray 
bands, triggering intensive observations in the radio, optical and 
X-ray (Swift, Chandra, INTEGRAL) bands (Villata at al. 2006, 
Giommi et al. 2006, Pian et al. 2006). Unfortunately no y-ray 
satellite was operating in the GeV domain at that time and no 
information was obtained in that band. 

During the summer of 2007, 3C454.3 was active again, 
reaching a level of the optical emission comparable to that 
of 2005. Several observations in the optical. X-ray and y-ray 
band were triggered (optical: KVA, optical-UV: Swift/UWOT, X- 
ray: Swift/XRT, GeV band: AGILE/GRID). The AGILE satellite 
(Tavani et al. 2008), still in its science verification phase, de- 
tected intense emission from 3C454.3 (Vercellone et al. 2008a). 

Triggered by these observations, the Major Atmospheric 
Gamma-ray Imaging Cherenkov (MAGIC) Telescope observed 
3C454.3 in July and August 2007. Another y-ray active 
phase was recorded by AGILE in November-December 2007 
(Vercellone et al. 2008b, Vercellone et al. 2009, in prep), which 
triggered further observations with MAGIC during that period. 
In the following (Section |2]i we describe the MAGIC observa- 
tions and the analysis procedure. In Sect. [3] we interpret the re- 
sults in the framework of the widely assumed Synchrotron Self 
+ External Compton (e.g. Sikora et al. 1994) model. 

2. MAGIC observations and data analysis 

MAGIC (Baixeras et al. 2004, Cortina et al. 2005) is a new gen- 
eration Imaging Atmospheric Cherenkov Telescope at La Palma, 
Canary Islands, Spain (28.3°N, 17.8°W, 2240 m asl). Thanks to 
its low energy trigger threshold of 60 GeV, MAGIC is well suited 
for multiwavelength observations together with the instruments 
operating in the GeV range. The parabolically-shaped reflector, 
with its total mirror area of 236 m^ allows MAGIC to sample a 
part of the Cherenkov light pool and focus it onto a multi-pixel 
camera, composed of 576 ultra-sensitive photomultipliers. The 
total field of view of the camera is 3.5° and the collection area is 
of the order of 10^ m^ at 200 GeV for a source close to zenith. 

The incident light pulses are converted into optical sig- 
nals, transmitted, via optical fibers and digitized by 2-GHz flash 
ADCs (Goebel et al. 2007). The primary particle energy and in- 
coming direction are reconstructed done by studying the inten- 
sity (and area) of the elliptical images and their orientation in the 
camera. In particular, the shape of the image allows to suppress 
hadron-induced showers and thus to reject the hadronic back- 
ground. 

In July and August 2007 observations were carried out in the 
ON/OFF mode, in which the source was observed on axis (for 
a total of 9.6 hours), while for background estimation, observa- 
tions (for a total of 7.3 hours) from a region of similar condi- 
tions in the sky from where no gamma rays are expected were 
used. Later, in November and December 2007, additional ob- 
servations were performed in the false-source tracking (wobble) 



mode (Fomin et al. 1994) in which the telescope was pointed al- 
ternatingly for 20 minutes to two opposite sky positions at 0.4° 
offset from the source (total 6.8 hours). This procedure allows 
to simultaneously determine the background and thus no extra 
OFF observations are needed. The zenith angle of all these ob- 
servations ranged from 12 to 30 degrees. The weather conditions 
in July were not as good as those of August, hence the event rate 
was lower. 

The analysis was performed using the standard MAGIC anal- 
ysis software (Bretz et al. 2005). After calibration and image 
cleaning based on a two-level tail cut (6 photoelectrons for im- 
age core and 3 photoelectrons for boundary pixels; see Fegan 
1997), the camera images were parameterized by the so-called 
Hillas image parameters (Hillas 1985). Two additional parame- 
ters, namely the time gradient along the main shower axis and 
the time spread of the shower pixels, were computed (Albert 
et al. 2008d). Hadronic background suppression was achieved 
using the Random Forest (RF) method (Breiman 2001, Albert 
et al. 2008c), in which for each event the so-called hadronness 
is computed, based on the Hillas and the time parameters. The 
"hadronness" parameter can be calculated for every event, and is 
a measure of the probability that the event is not y like. The RF 
method was also used for the energy estimation. Crab Nebula 
data from the same periods and zenith angle distributions were 
studied using the same analysis chain to check the validity of the 
results. 

Since there was no significant signal found, upper limits 
(95% CL) were calculated (Rolke et al. 2005) taking into ac- 
count a 30% systematic error in energy determination and effec- 
tive area calculation (see Albert et al. 2008b). Table[T]shows the 
results for the July-August observations, whereas upper limits 
for the November-December observation are given in Tabled 
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U.L. July 17-20 2007 




[GeV] 


NoN 


NoFF 


Sign, (cr) 


C.U. 


[erg cm - s '] 


83 


54188 


54705 


-1.56 


0.04 


0.78 X 10-" 


186 


976 


965 


0.25 


0.05 


0.62 X 10-" 


476 


62 


52.3 


0.91 


0.03 


0.169 X 10-" 


<£> 




U.L. August 9-22 2007 


[GeV] 


NoN 


NoFF 


Sign, (cr) 


C.U. 


[erg cm - s ] 


128 


5453 


5539.7 


-0.82 


0.14 


2.0 X 10-" 


186 


3892 


3885.1 


0.078 


0.03 


0.3 X 10-" 


476 


202 


220.8 


-0.91 


0.01 


0.09 X 10-" 



Table 1. Derived upper limits on flux for the July and the August 
2007 data. The columns represent respectively: the average true 
energy, the number of ON source events , number of background 
(OFF) events, the significance, the flux upper limit in Crab Units 
(C.U.) and in absolute flux units of [erg cm-^ s-']. 



3. Discussion 

The Spectral Energy Distribution (SED) of 3C454.3 around the 
epoch of the July and November 2007 MAGIC observations, as- 
sembled with the available data, is shown in Figure [1] For July 
(left panel) we show the nearly simultaneous data in the opti- 
cal (KVA, July 24), optical-UV (Swift/UYOT, July 26), X-ray 
(SwiftpCKT, July 26) and y-ray (AGILE/GRID) band (average of 
July 24-30). For November (right panel) the data in the optical. 
X-ray (Swift/XRT and INTEGRAL) and y-ray (AGILE/GRID) 
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<£> 
[GeV] 


NoN 


U.L.I 

NoFF 


Vov.27,30 & 

Sign, (o-) 


Dec. 1 

C.U. 


2007 

[erg cm"^ s 


113 


39900 


39920 


-0.07 


0.3 


4.6 X 10^" 


235 


385 


367 


0.66 


0.09 


0.9 X 10^" 



Table 2. Derived upper limits on flux for the November and the 
December 2007 data. The columns represent respectively: the 
average true energy, the number of ON source events, number of 
background (OFF) events, the significance, the flux upper limit 
in Crab Units (C.U.) and in absolute flux units of [erg cm"^ s"']. 



data, averaged over the entire period of the AGILE observa- 
tions (Nov. I I-Dec. 1) are shown. For comparison we also show 
(open circles) historical data. AGILE/GRID spectra, both for the 
July and November observations, have been recently published 
in Vercellone et al. (2008b). 

In the same figure, upper limits from MAGIC observations 
(18-21 July and 27, 28, 30 November) are shown as triangles 
(observed: empty; EBL-deabsorbed: filled) (see Table For 
the EBL deabsorption we used the LowSFR model of Kneiske 
et al. (2004) which predicts a low level of the EBL close to 
what is presently inferred from observations, both directly (e.g. 
Franceschini et al. 2008) and indirectly (Aharonian et al. 2006b, 
Mazin & Raue 2007, Albert et al. 2008a). 

Fig. [H shows that the (absorption-corrected) MAGIC upper 
limit at ~100 GeV is inconsistent with the extrapolation of the 
hard y-ray (100 MeV-10 GeV) spectrum. Therefore, the data in- 
dicate there is a break (or a cutoff) of the emission between the 
GeV and the 100 GeV band. As discussed below this is consis- 
tent with the expectations from the simplest leptonic model. 

Emission from blazars is dominated by the non-thermal con- 
tinuum emitted by a relativistic jet closely aligned towards the 
observer. The SED of FSRQs is widely interpreted in terms of 
synchrotron and inverse Compton emission from high-energy 
electrons. The latter component is probably dominated by the 
scattering of the external photons (originating in the disk and/or 
in the broad line region [BLR], Sikora et al. 1994), although the 
synchrotron self-Compton emission (Maraschi et al. 1992) and 
the inverse Compton scattering of the direct radiation from the 
accretion disk (e.g. Dermer & Schlickeiser 1993) can signifi- 
cantly contribute in the X-ray band. The SED of 3C454.3, in- 
cluding optical. X-rays and GeV measurements around the end 
of July, has been already discussed and modelled by Ghisellini et 
al. (2007). However, the model discussed in that work assumes 
that the y-ray spectrum was similar to the average EGRET spec- 
trum, with a soft slope. The spectrum of AGILE (Vercellone et 
al. 2008b), instead, is rather hard (photon index F ^ 1 .7 in the 
100 MeV-1 GeV band), both for July and November, suggest- 
ing a peak of the high-energy component at frequencies above 
~ 10^^ Hz. The SEDs of November has been already discussed 
and modeled in Vercellone et al. (2008b). 

To reproduce the multifrequency data we use the emission 
model fully described in Maraschi & Tavecchio (2003). Given 
the focus on the VHE emission, we also consider the absorp- 
tion of y-ray photons through pair production within the BLR. 
Moreover, the external radiation field (assumed to be isotropic 
in the frame of the black hole), usually approximated by a black 
body emission peaking in the UV region, has been calculated 
using the photoionization code CLOUDY (Ferland et al. 1998). 
Details on the emission model can be found in Maraschi & 
Tavecchio (2003), while the description of the calculation of 
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Table 3. Parameters used in the emission model to reproduce the 
SEDs of 3C454.3 in July (first row) and November (second row) 
2007. F: bulk Lorentz factor of the emission region; B: magnetic 
field, in units of [G]; K: density of the relativistic electrons, in 
units of [cm^^]; n\ and low and high energy slope of the 
electron energy distribution; jmrn, Jb and ymax^ minimum energy, 
break energy and maximum energy of the electron energy distri- 
bution; R: radius of the emission region, in units of [10'^ cm]. 
See text for more details. 



the external radiation field is reported in Tavecchio & Ghisellini 
(2008). 

We assume that the emission is produced within a spheri- 
cal region of radius R, in motion with bulk Lorentz factor F. 
We assume that the corresponding relativistic Doppler factor is 
5 - I IT. The tangled magnetic field has an intensity B. The 
emitting particles, with density K, follow a (steady state) broken- 
power law energy distribution extending from y^ia to 7max, with 
indices n\ and «2 below and above the break at jh- This purely 
phenomenological distribution has been assumed to reproduce 
the observed shape of the blazar SEDs, without any specific as- 
sumption on the acceleration/cooling mechanism acting on the 
particles. With this choice we are allowed to assume extreme 
low-energy slopes («i < 2) such as those required for 3C454.3, 
which cannot be obtained under standard conditions. It is con- 
ceivable that, at least in these cases, the electron distribution 
derives from two (continuously operating) diff'erent acceleration 
mechanisms (see e.g. Sikora et al. 2002). We also neglect the 
effects related to the cooling of particles in the Klein-Nishina 
regime, discussed by Moderski et al. (2005). We note, how- 
ever, that these effects should produce a bump in the optical- 
UV synchrotron emission which is not apparent in the available 
data, though the poor coverage does not allow a firm conclu- 
sion. We model the external radiation field assuming that the 
disk emission (dotted line in figure), with a total luminosity of 
idisk = 5 X 10*** erg/s, is reprocessed by clouds of the BLR, a 
sphere with radius 3 x 10'^ cm (we assume that clouds are char- 
acterized by standard values of the density hblr =10" cm"^ and 
hydrogen column density, A^h = 10^^ cm"^). For simplicity, we 
assume that the distance of the emission region of the jet from 
the central black hole is smaller than the radius of the BLR, but 
large enough to neglect the direct disk emission, coming from 
behind the jet (e.g. Dermer & Schlickeiser 1993; Vercellone et 
al. 2008b consider in the model also this component). To repro- 
duce the shape of the high-energy component, we assume that 
in the X-ray band the emission is dominated by SSC emission, 
while EC radiation accounts for the GeV peak. 

In our model, the rapid decrease of the flux above few tens of 
GeV would be related to two effects: i) the decrease of the scat- 
tering cross section and ii) the absorption of the produced y-rays 
through pair production. The energy above which the KN effects 
become important can be roughly expressed as: /skn - 22.5v^ 
GeV, where v,, 15 is the frequency of the external photons (in 
units of 10'^ Hz). The emission including only the KN effects, 
neglecting the absorption, is shown by the long dashed-dotted 
line. The frequency above which the absorption of y-rays be- 
come effective can be roughly expressed as: Eyy ^ fiOv^J^ GeV, 
where v,,i5 is the frequency of the target photons (in units of 
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Fig. 1. SED of 3C454.3 assembled with multifrequency information available for the period close to the MAGIC observation at the 
end of July 2007 (left panel; optical: KVA, optical-UV: Swift/VWOT, X-ray: Swift/XRT, GeV band: AGILE/GRID) and November 
2007 (right panel: optical-UV: Swift/UVOT, X-ray: Swift/XRT and INTEGRAL, GeV band: AGILE/GRID). Triangles report the 
observed (empty) and the deabsorbed (filled) upper limits of MAGIC in three different bands. For comparison we also report (open 
circles) historical data (Kuhretal. 1981, NED, Gearetal. 1994, Stevens etal 1994, Impey & Neugebauer 1988, Smith et al. 1988 for 
radio and optical; Tavecchio et al. 2007 for X-rays from Chandra). The open circle and the bow-tie in the MeV-GeV region indicate 
the average EGRET spectrum (Hartman et al. 1999). Solid line reports the results of the modelling with the synchrotron-inverse 
Compton model (see text for details and model parameters). We also report the single emission components: synchrotron (dashed), 
SSC (dotted-dashed) and EC (long dashed). The dotted line shows the emission of the accretion disk. 



10'^ Hz). Therefore, as shown by the solid line in Fig [T] (calcu- 
lated including both effects, the internal absorption treated as in 
Tavecchio & Mazin 2008), the expected emission above 20-30 
GeV is rather small, consistent with the observed upper limits. 
Note that, although the limit set by KN effects is a characteris- 
tic feature of leptonic models, absorption of y-rays by soft pho- 
tons can also be relevant for hadronic models (e.g. Reimer 2007). 
More detailed discussions on the effects of absorption, including 
the possible hardening of the observed spectrum above ~ 100 
GeV, can be found in, e.g., Aharonian et al. (2008b), Sitarek & 
Bednarek (2008), Tavecchio & Mazin (2008). 

Summarizing, we have shown that the upper limits in the 
VHE band for the 3C454.3 obtained with the MAGIC telescope 
are consistent with the expectations of the leptonic models for 
FSRQs, predicting a sharp decrease of the flux above few tens 
of GeV, due to the internal absorption of y-rays and the de- 
creased efficiency of the inverse Compton emission at high en- 
ergy. Therefore even upper limits, particularly if accompanied by 
simultaneous observations in the MeV-GeV band, can be useful 
to test current emission models for FSRQs. Stronger constraints 
will definitely be obtained with future multifrequency cam- 
paigns already planned with the Fermi Gamma-Ray Telescope. 
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